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ABSTRACT 


In  1998,  the  Naval  Postgraduate  School  (NPS)  Center  for  Interdisciplinary  Remotely 
Piloted  Aircraft  Studies  (CIRPAS)  and  the  NPS  Department  of  Electrical  and  Computer 
Engineering  (ECE  Dept.)  collaborated  on  the  acquisition  of  a  mobile  radar,  the  AN/MPQ- 
64  (Sentinel).  This  is  a  modern  X-band,  pulse  Doppler  radar  used  by  the  Army  for 
forward  area  air  defense.  An  SBIR  project  funded  by  the  Office  of  Naval  Research 
resulted  in  a  contract  to  ProSensing,  Inc.,  Amherst,  MA  to  retrofit  this  radar  with  a 
weather  processor.  The  intent  is  to  provide  a  military  capability  to  assess  battlespace 
weather  in  real  time.  This  capability  will  be  developed  using  the  NPS  radar  as  a  testbed. 
When  the  weather  capability  has  been  fully  implemented  on  the  NPS  radar,  it  will  also 
have  application  as  a  scientific  instrument  for  severe  storm  research  because  of  its 
mobility  and  cutting-edge  capabilities.  The  modified  radar  has  been  designated  the 
MWR-05XP  (Mobile  Weather  Radar,  2005,  X-band,  Phased  Array). 

This  report  presents  the  results  of  an  analysis  of  numerous  aspects  of  the  radar’s 
performance  as  a  weather  sensor.  The  ability  of  the  MWR-05XP  radar  to  detect  rain 
with  different  levels  of  reflectivity,  the  ability  of  the  radar  to  detect  clear  air  turbulence 
with  different  structure  parameters  and  echoes  from  birds  and  insects  has  been  examined. 
Curves  are  provided  to  determine  correlation  and  decorrelation  times  for  weather  signals 
with  varying  rms  velocity  spread.  MWR-05XP  post  detection  integration  improvement 
has  been  computed  for  a  Rayleigh  weather  target  with  varying  pulse-to-pulse  correlation 
and  curves  are  provided  to  determine  the  improvement  in  terms  of  the  number  of  pulses 
integrated.  Lastly,  scan  strategy  is  discussed  with  emphasis  on  obtaining  weather  signal 
parameter  estimates  with  small  variance  while  at  the  same  time  achieving  a  rapid 
volumetric  update  rate.  Electronic  scan  is  the  feature  that  enables  the  radar  to  achieve  a 
rapid  volumetric  update  rate.  Data  quality  control  is  also  discussed. 


(Back  to  Table  of  Contents) 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


IV 


EXECUTIVE  SUMMARY 


Theoretical  results  provided  in  this  report  predict  the  performance  of  an  electronically 
scanned  X-band  radar  equipped  with  a  weather  signal  processor.  The  radar  is  the  MWR- 
05XP  (Mobile  Weather  Radar,  2005,  X-band,  Phased  array).  The  single  pulse  average 
power  return  from  precipitation,  birds,  insects,  flocks  of  birds,  swarms  of  insects  and 
clear  air  (turbulence  &  wind  shear)  has  been  investigated.  A  postdetection  integration 
gain  curve  keyed  to  the  MWR-05XP  is  provided  along  with  a  curve  of  single  pulse 
fluctuation  loss.  Integration  gain  curves  for  weather  signal  samples  with  correlation 
values  from  0.40  to  0.99  are  provided.  These  curves  quantify  integration  gain  as  a 
function  of  the  number  of  pulses  integrated  taking  signal  fluctuation  and  pulse-to-pulse 
correlation  into  account.  Lastly,  an  optimum  sampling  strategy  is  proposed.  The 
proposed  strategy  maximizes  the  volumetric  update  rate  while  simultaneously 
minimizing  the  variance  of  signal  parameter  estimates. 

The  MWR-05XP  is  designed  for  a  maximum  operational  range  of  approximately  75  km. 
A  PRF  of  2  kHZ  is  required  to  achieve  an  unambiguous  range  of  75  km.  For  the 
operating  frequency  range  of  the  MWR-05XP,  this  PRF  results  in  an  unambiguous 
velocity  of  about  30  m/s.  Increasing  PRF  to  the  system  limit,  10  kHz,  increases  the 
unambiguous  velocity  to  about  150  m/s  while  unambiguous  range  drops  to  15  km.  The 
latter  limits  are  well  suited  for  use  of  the  MWR-05XP  as  a  weather  research  instrument 
where  the  radar  would  probably  be  positioned  to  observe  severe  weather  at  short  range. 
The  former  limits  would  be  suitable  for  integration  of  an  ordinary  weather  observation 
capability  into  the  normal  military  operational  mode  of  the  radar. 

Single  pulse  received  signal  power  for  precipitation  has  been  computed  for  reflectivity 
levels  between  5  dBZ  and  75  dBZ  and  ranges  from  5  to  45  km.  These  computations 
show  average  echo  signal  power  is  generally  in  the  range  between  the  system  noise  floor 
and  receiver  saturation  level.  The  reflectivity  range  5-50  dBZ  corresponds  to  rainfall 
rates  of  approximately  0.06-37.7  mm/hr.  Larger  reflectivities  are  likely  due  to  melting 
ice.  The  results  show  observation  of  severe  weather  with  high  reflectivity  at  short  range 
is  likely  to  cause  receiver  saturation.  Even  for  reflectivities  with  average  power  return 
below  saturation,  individual  pulses  may  drive  the  receiver  into  saturation  as  they  can 
typically  exceed  the  average  by  as  much  as  10  dB. 

In  the  absence  of  precipitation,  return  can  be  observed  from  birds,  insects  and  clear  air. 
Computations  show  the  X-band,  single  pulse  echo  signal  power  from  individual  birds  and 
large  insects  such  as  moths  will  he  above  the  system  noise  floor  if  range  is  less  than 
about  10  km.  Flocks  of  birds  and  insect  swarms  produce  a  much  greater  return  with 
average  power  very  similar  to  that  for  echoes  from  precipitation.  In  fact,  research  has 
shown  that  “clear  air”  return  is  frequently  associated  with  insects  and  has  a  diurnal 
pattern.  Birds,  although  they  are  targets  with  larger  cross-section,  are  not  as  numerous  so 
their  RCS  density  is  generally  less  than  that  of  insects. 
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True  clear  air  return  results  from  variations  in  refractive  index  due  to  turbulence  and  wind 
shear.  Turbulence  is  the  predominant  mechanism.  The  return  from  clear  air  turbulence  is 
weak,  however,  and  computations  show  single  pulse  echo  signal  power  is  typically  more 
than  20  dB  below  the  MWR-05XP  system  noise  floor  even  for  strong  turbulence.  Thus, 
it  will  not  normally  be  a  factor. 

If  clear  air  turbulence  is  of  interest,  the  signal  can  be  processed  if  a  sufficiently  long 
observation  time  is  used.  Computations  show  strong  turbulence  at  a  range  of  10  km 
produces  a  single  pulse  echo  signal  with  power  P.  «  -134  dbm  in  X-band.  An 

integration  gain  of  2 1  dB  would  bring  this  signal  to  the  system  noise  floor  and  further 
integration  would  produce  SNR  >  0  dB  so  the  signal  could  be  detected.  Several  thousand 
pulses  would  be  required  and  this  would  consume  a  time  on  the  order  of  one  second  for 
each  beam  position.  A  time  on  the  order  of  minutes  would  be  required  to  scan  a  frame  of 
reasonable  size.  For  reference,  the  S-band  WSR-88D  NEXRAD  radar  displays  clear  air 
reflectivity  between  -28  dBZ  and  +28  dBZ  using  a  10  minute  observation  time  for  its 
volume  scan.  The  volume  scan  covers  360  degrees  in  azimuth  and  0.5  to  4.5  degrees  in 
five  elevation  beam  steps.  The  lowest  two  elevation  angles  are  scanned  twice  using 
different  pulse  widths  and  PRFs  for  reflectivity  and  velocity. 

Pulse  integration  can  be  used  as  desired  to  improve  signal-to-noise  ratio,  SNR.  The 
signal  from  precipitation  is  noncoherent  on  a  relatively  short  time  scale  so  noncoherent 
integration  must  be  used.  Albersheim’s  equations  have  been  used  to  model  the 
noncoherent  integration  gain,  Gnc ,  for  an  arbitrary  number  of  integrated  pulses.  The 

result  also  depends  on  detection  and  false  alarm  probabilities.  Pd  =  0.95  and  Tfa  =  1  hour 

have  been  chosen  as  reasonable  for  the  MWR-05XP.  For  this  choice  of  parameters, 
integration  of  500  pulses  produces  a  noncoherent  integration  gain,  Gnc  =  19  dB  . 

Precipitation  is  a  Rayleigh  target,  the  observed  radar  cross-section  is  exponentially 
distributed  and  thus  echo  signal  power  varies  with  time.  The  echo  signal  power  for  any 
particular  pulse  might  be  significantly  less  than  the  average  echo  signal  power.  There  is  a 
resulting  fluctuation  loss  that  impacts  single  pulse  detection  probability.  This  fluctuation 
loss  has  been  modeled  using  Shnidman’s  equations.  For  a  single  pulse  detection 
probability  Pd  =  0.95  the  fluctuation  loss  is  12  dB.  For  a  single  pulse  detection 

probability  Pd  =  0.99  the  fluctuation  loss  is  17.9  dB.  This  loss  must  frequently  be 
overcome  by  using  integration. 

Work  by  Kanter  showed  there  will  be  a  reduction  in  fluctuation  loss  if  pulses  are 
integrated.  The  reduction  in  fluctuation  loss  when  integrating  n  pulses  can  be  modeled 
using  an  equivalent  number  of  pulses,  ne,  that  depends  on  correlation  coefficient. 
Computations  have  been  carried  out  and  results  show  that  when  correlation  coefficient 
p  >  0.4  ,  the  effective  number  of  pulses  is  less  than  the  actual  number.  As  p  — »  1  (pulses 
completely  correlated)  the  effective  number  of  pulses,  ne  — »  1 . 
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Integration  improvement  can  be  expressed  in  terms  of  integration  gain  and  fluctuation 
loss  for  pulses  with  arbitrary  correlation  coefficient.  Computations  have  been  carried  out 
and  curves  are  provided  that  quantify  the  integration  improvement  for  correlation 
coefficients  ranging  from  p  =  0.4  to  p  =  0.99  .  These  curves  show  that  if  a  sufficient 
number  of  pulses  are  integrated,  the  effect  of  weather  target  cross-section  fluctuations  is 
masked.  The  integration  improvement  approaches  that  for  a  non-fluctuating  target  with 
RCS  <J^=E{a}. 

The  last  issue  considered  in  this  work  was  optimum  sampling  strategy.  Weather  radars 
typically  produce  estimates  of  the  first  three  Doppler  moments  of  the  weather  signal; 
reflectivity,  average  radial  velocity  and  radial  velocity  spread.  The  estimates  have  a 
variance  that  depends  on  the  number  of  samples  averaged  and  their  correlation 
coefficient,  p.  For  a  given  number  of  samples,  p  =  0  yields  the  estimate  with  minimum 
variance.  An  optimum  sampling  strategy  is  defined  here  as  one  that  maximizes  the 
volumetric  update  rate  while  simultaneously  minimizing  the  variance  of  the  weather 
signal  parameter  estimates. 

Reflectivity  is  estimated  from  the  average  power  return.  Average  radial  velocity  is 
estimated  using  the  average  of  phase  shifts  computed  between  pulse  pairs  (the  pulse  pair 
algorithm).  Pairs  of  pulses  must  be  coherent  (highly  correlated)  and  therefore  closely 
spaced  in  time  to  unambiguously  compute  velocity  from  the  inter  pulse  phase  shift. 
However,  pulse  pairs  should  be  widely  spaced  in  time  if  successive  power  and  velocity 
estimates  are  to  be  independent  (uncorrelated). 

The  Doppler  spectrum  of  a  weather  signal  is  assumed  Gaussian  with  spectral  spread 
directly  tied  to  weather  target  rms  velocity  spread,  av.  The  correlation  coefficient  of  the 
signal  is  found  from  the  Fourier  transform  of  the  power  spectral  density.  Therefore,  the 
correlation  coefficient  of  the  signal  is  also  Gaussian.  The  correlation  coefficient  for  an 
X-band  weather  signal  has  been  computed  for  a  typical  range  of  nns  velocity  spreads, 

0. 1  <  cr ,  <  10  m/s  .  The  correlation  coefficient  has  been  used  to  determine  the  time  the  X- 

band  signal  remains  correlated  as  well  as  the  time  required  for  the  signal  to  decorrelate. 
These  times  detennine  the  maximum  acceptable  spacing  between  pulses  in  a  pair  and  the 
time  one  must  wait  to  obtain  another  independent  sample  pair  from  the  same  resolution 
cell. 

An  optimum  sampling  strategy  follows  from  the  above  considerations.  A  pair  of  pulses 
should  be  transmitted  to  collect  a  single  sample  pair  for  each  resolution  cell  of  interest  in 
a  beam  position.  The  process  is  repeated  for  each  of  the  beam  positions  in  a  frame.  After 
sample  pairs  have  been  collected  for  all  beam  positions  in  a  frame,  the  elapsed  time 
should  exceed  the  weather  signal  decorrelation  time  and  the  process  can  be  repeated  to 
collect  an  independent  set  of  sample  pairs. 

Implementation  of  the  optimum  sampling  strategy  requires  electronic  antenna  beam  scan. 
Sample  time  for  a  beam  position  is  detennined  by  the  spacing  between  pulses  in  a  pair 
and  the  time  required  to  electronically  switch  beam  positions.  The  time  required  for  the 
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MWR-05XP  to  switch  beam  positions  is  nominally  100  ps.  At  the  maximum  PRF  of  10 
kHz,  the  correlation  of  adjacent  pulses  is  high  and  the  time  between  pulses  is  100  ps,  the 
same  as  the  nominal  phase  shifter  switching  time.  The  total  time  required  to  collect  a 
sample  pair  and  switch  beam  positions  is  thus  300  ps.  Sample  pairs  could  be  collected 
from  200  beam  positions,  for  example,  in  60  ms.  After  60  ms  has  elapsed,  the  weather 
signal  decorrelates  and  another  scan  of  the  frame  yields  independent  sample  pairs. 
Unambiguous  range  for  this  scenario  would  be  15  km. 

It  is  shown  that  the  sampling  strategy  described  above  allows  for  simple  control  of  data 
quality.  The  desired  variance  of  a  parameter  estimate  can  be  specified  and  since  samples 
are  independent,  the  number  of  samples  required  can  be  easily  detennined.  Collecting 
the  computed  number  of  samples  with  sufficiently  high  SNR  then  results  in  weather 
signal  parameter  estimates  with  known  variance.  In  the  case  of  reflectivity,  however,  the 
radar  must  be  accurately  calibrated  so  echo  signal  power  at  the  receiver  input  can  be 
determined  from  the  observable  I  and  Q  channel  output  voltages.  Each  dB  of  calibration 
error  will  bias  the  reflectivity  estimate  by  one  dBZ. 
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I.  INTRODUCTION 


A.  BACKGROUND 

The  Naval  Postgraduate  School  (NPS)  Department  of  Electrical  and  Computer 
Engineering  (http://www.nps.navy.mil/cce/),  the  NPS  CIRPAS  Research  Center 
(http://web.nns.navy.mil/~cirpas/)  and  ProSensing,  Inc.  (http://www.prosensing.com/) 
have  collaborated  on  the  development  of  a  mobile,  phased  array,  pulse  Doppler 
weather  radar.  The  system  was  created  from  an  AN/MPQ-64  mobile  radar  obtained 
from  the  Army  (Figure  1).  A  Weather  Radar  Processor  (WRP)  was  developed  by 
ProSensing,  Inc.  (http://www.prosensing.com/)  under  a  Phase  II  SBIR  project  (Topic 
N01-035)  funded  by  the  Office  of  Naval  Research  (ONR).  The  weather  processor 
was  added  to  the  radar  and  the  radar  was  mounted  on  a  flatbed  truck  (Figures  2,  3). 

NPS  has  a  Cooperative  Research  and  Development  Agreement  with  ProSensing,  Inc. 
and  has  been  collaborating  with  ProSensing  on  research  related  to  the  future 
operational  use  of  the  radar  as  a  weather  instrument.  In  keeping  with  the 
nomenclature  applied  to  the  National  Weather  Service  NEXRAD  radars  (WSR-88D, 
Weather  Surveillance  Radar,  1988),  the  author  has  designated  the  radar  described  in 
this  report  the  MWR-05XP  (Mobile  Weather  Radar,  2005,  X-band,  Phased  Array). 
The  radar  can  be  electronically  scanned  and  it  is  this  feature  that  gives  it  a  unique 
capability  as  a  weather  research  instrument.  That  capability  is  rapid  volumetric 
update  of  weather  signal  parameter  estimates  while  simultaneously  achieving  small 
variance  of  the  estimates.  /Back  to  Table  of  Contents) 

B.  OBJECTIVE 

The  objective  of  the  work  described  herein  was  to  investigate  a  number  of  issues 
related  to  the  operation  of  the  MWR-05XP  as  a  weather  sensor.  These  included  range 
and  velocity  ambiguity,  sensitivity,  return  from  precipitation,  clear  air  return,  bird  and 
insect  return,  fluctuating  target  integration  improvement  and  scan  strategy.  This 
report  should  serve  as  a  useful  reference  on  the  performance  of  the  MWR-05XP. 
(Back  to  Table  of  Contents) 

C.  APPROACH 

The  approach  throughout  is  analytical  and  builds  on  existing  theory  to  develop  results 
specific  to  the  operation  of  the  MWR-05XP  as  a  weather  sensor.  Basic  radar  system 
measurements  were  made  to  support  the  analysis;  Appendix  2.  (Back  to  Table  of 
Contents) 

D.  RELATED  WORK 

There  is  a  significant  body  of  literature  relating  to  the  use  of  radar  for  weather 
measurements.  Much  of  this  relates  to  the  WSR-88D  NEXRAD  weather  radar. 
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Figure  1.  MWR-05XP  (Mobile  Weather  Radar-2005,  X-band,  Phased  Array).  Radar 
shown  in  its  original  military  operational  configuration,  mounted  on  a  trailer.  (Back) 


Figure  2.  Drawing  showing  MWR-05XP  radar  mounted  on  a  flatbed  truck  with  400 
Hz  diesel  generator  and  operators’  shelter.  Antenna  is  in  stowed  position.  Radar  is 
mounted  on  a  platform  used  to  raise  the  radar  above  the  operator’s  shelter  when  in 
use.  Courtesy  ProSensing,  Inc.  (Back) 
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Figure  3.  Photograph  of  MWR-05XP  mobile  weather  radar  shown  with  radar  platform 
raised  for  operation.  The  main  400  Hz  diesel  generator  on  the  rear  of  the  truck  powers 
the  radar.  A  120  V,  60  Hz  generator  on  the  driver’s  side  below  the  truck  bed  provides 
power  for  computers  and  equipment  in  the  operator’s  shelter.  (Back) 

There  are  approximately  160  NEXRAD  radars  that  provide  coverage  of  the  continental 
US,  Alaska,  Hawaii,  Guam  and  Puerto  Rico.  These  radars  can  be  accessed  on  the 
internet  using  the  URL  (http://weather.noaa.gov/radar/national.html).  Doppler  radar 
theory  and  meteorology  are  described  in  Ref.  [1].  A  rigorous  treatment  of  Doppler  radar 
and  weather  observations  can  be  found  in  the  text  by  Doviak  and  Zrnic,  Ref.  [2].  Other 
texts  covering  weather  radar  theory  include  Sauvageot,  Ref.  [3],  and  Rinehart,  Ref.  [4]. 
Skolnik,  Ref.  [5],  also  discusses  radar  return  from  weather  as  well  as  land  clutter,  sea 
clutter,  bird  and  insect  return.  Surface  return  is  clutter  for  a  weather  radar  but  that 
problem  is  not  treated  here.  Martin  [6]  has  investigated  return  from  birds  and  insects. 
Others  [7]  -  [12]  have  published  work  related  to  the  detection  of  radar  targets  and  in 
particular,  the  detection  of  signals  from  fluctuating  targets  such  a  precipitation.  These 
theoretical  results  have  been  used  to  predict  the  performance  of  the  MWR-05XP  weather 
radar  and  to  determine  temporal  and  spatial  sampling  strategies  that  fully  exploit  the 
electronic  scan  capability  of  the  radar.  (Back  to  Table  of  Contents) 
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II.  ANALYSIS  OF  MWR-05XP  PERFORMANCE 


A.  WEATHER  RADAR  EQUATION 


The  weather  radar  equation  can  be  found  in  various  forms  in  a  number  of  sources  (Refs. 
[1]  -  [5])  and  will  not  be  derived  here.  The  average  single  pulse  return  echo  signal  power 
from  precipitation  is  given  by 

'■  1024ln(2)tf2L,  /  (  ’ 

where 

Pr  =  average  received  echo  signal  power  in  Watts 

Pt  =  transmit  power  in  Watts 

G,  =  transmit/receive  antenna  gain 

X  =wavelength  in  meters 

6b  =  antenna  3  dB  elevation  beamwidth 

(pB  =  antenna  3  dB  azimuth  beamwidth 
c  =  3x  10s  m/s 

|AT|~  =|(fr  -l)/(fr  +2)|2  «  0.93  (water) 

R  =  range  in  meters 
Ls  =  system  loss 
Z=  reflectivity  in  mmV. 

The  nature  of  a  weather  target  is  such  that  pulses  are  normally  integrated  noncoherently 
(post-detection).  If  pulses  are  integrated,  then  an  integration  gain  factor,  Gnc(n )  and 
fluctuation  loss  factor  must  be  added  [5],  This  yields 

-  _  (e,<p. A2)(cr)  \K\2  „ 

r  .  .  1/  T  4  ^ 

1 024  In  (2)  R2Ls(Lfyne  A 

where 

Lf=  single  pulse  fluctuation  loss  factor 
n  =  number  of  pulses  integrated 
ne(n,p)  =  min  [n,l  +  (n-l)\n(l/p)] 
p  =  correlation  coefficient . 

Equation  (2)  is  the  basis  for  computation  of  the  basic  weather  detection  capabilities  of  the 
MWR-05XP. 
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It  is  useful  introduce  signal-to-noise  ratio  in  Eq.  (2).  This  is  accomplished  simply  by 
dividing  the  left  side  by  noise  power  N,  and  the  right  side  by  noise  power  kT  B  to  obtain 


lp/\  PtfG«(n)X-(9ji>.ln2)(cT)  ejxf 


V/Af 


10241n(2)^(iJ;S„)4(i/) 


V 


X' 


(3) 


where 


A:  =  1.38x10  23  W/Hz-°K 
Ts  =  system  noise  temperature  in  °K 
Bn  =  equivalent  noise  bandwidth 
pJn  =  mean  SNR. 


(Back  to  Table  of  Contents) 


B.  UNAMBIGUOUS  RANGE  AND  VELOCITY 

The  relations  for  unambiguous  range  and  unambiguous  velocity  are  given  by  [1]  -  [5] 

Ru  = 

and 

v,  = 


c 

(3) 

2fp 

V, 

4 

(4) 

where 


R„  =  unambiguous  range  in  meters 
vu  =  unambiguous  velocity  in  m/s 
fp  =  pulse  repetition  frequency  in  Hz 
X  =  wavelength  in  meters 
c  =  3xl0N  m/s . 

Thus,  the  selection  of  PRF  is  a  tradeoff  when  both  target  range  and  velocity  are  of 
interest.  This  is  illustrated  in  Figure  4.  (Back  to  Table  of  Contents) 
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MWR-05XP  Unambiguous  Range  &  Velocity 
fi=9.4  GHz,  f2=10  GHz 

Prof.  Jeffrey  B.  Knorr 

Naval  Postgraduate  School  Monterey,  CA  93943 
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Figure  4.  Unambiguous  range  and  unambiguous  velocity  vs.  pulse  repetition  frequency 
for  the  MWR-05XP  radar.  (Back) 


It  will  be  shown  later  that  a  PRF  above  2500  Hz  is  generally  desired  for  MWR-05XP 
velocity  estimation.  Figure  4  shows  that  unambiguous  velocity  is  about  50  knots  and 
unambiguous  range  is  about  50  km  for  a  PRF  of  3  kHz.  For  a  PRF  of  10  kHz, 
unambiguous  velocity  can  be  increased  to  about  150  knots  but  unambiguous  range  drops 
to  15  km.  (Back  to  Table  of  Contents) 


C.  PRECIPITATION  SENSITIVITY 


The  sensitivity  of  a  weather  radar  is  often  shown  by  computing  and  displaying  the 
reflectivity  level  in  dBZ  that  the  radar  can  detect  at  a  specified  range  using  a  single  pulse 
with  SNR  =  0  dB.  Figure  5  depicts  the  sensitivity  of  the  MWR-05XP  weather  radar.  The 
figure  shows  that  a  reflectivity  of -5  dBZ  can  be  detected  at  a  range  of  10  km  while  for  a 
range  of  75  km,  the  detectable  level  of  reflectivity  rises  to  about  12.5  dBZ.  This  is  quite 
reasonable  as  a  scale  of  5  -  75  dBZ  is  normally  used  to  display  the  weather  return  when 
precipitation  is  present.  A  return  of  12.5  dBZ  corresponds  to  a  very  light  rainfall  rate, 
approximately  0.13  mm/hr. 

It  is  perhaps  of  interest  to  compare  the  sensitivity  of  the  MWR-05  XP  weather  radar  with 
that  of  the  WSR-88D  NEXRAD  radar.  The  latter  radar  operates  at  S-band.  The  longer 
wavelength  is  a  relative  disadvantage  for  scattering  in  the  Rayleigh  region  but  this 
disadvantage  is  overcome  through  use  of  a  higher  gain  antenna  (45  dB)  and  higher 
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Figure  5.  Single  pulse  reflectivity  sensitivity  in  dBZ  vs.  range  in  km  for  the  MWR-05XP 
radar.  SNR  =  0  dB.  ZdB/  =  -27.9  +  20  log  Rkm .  (Back) 


transmit  power  (750  kW).  The  WSR-88D  also  has  a  longer  pulse  width  than  that  of  the 
MWR-05XP.  Figure  6  shows  the  sensitivities  of  both  the  MWR-05-XP  and  the  WSR- 
88D.  It  can  be  seen  that  the  the  higher  ERP  and  lower  noise  floor  of  the  WSR-88D  radar 
more  than  offset  the  disadvantage  of  longer  wavelength,  resulting  in  a  net  sensitivity 
advantage  of  about  19  dB  for  the  WSR-88D.  Although  this  comparison  is  useful  because 
the  WSR-88D  is  a  well-known  system,  it  should  be  remembered  that  the  applications  and 
antenna  scans  of  the  two  radars  are  quite  different.  The  MWR-05XP  enjoys  its  own 
advantages,  mobility  and  electronic  scan.  (Back  to  Table  of  Contents) 
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MWR-05XP  vs.  WSR-88D  Single  Pulse  Sensitivity  for  SNR  =  0  dB 

Prof.  Jeffrey  B.  Knorr 

Naval  Postgraduate  School  Monterey,  CA  93943 
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Figure  6.  Single  pulse  reflectivity  sensitivity  comparison  of  the  MWR-05XP  and  WSR- 
88D  weather  radars.  ZfBZ'  =  -27.9  + 1 0  log  Rkm  .  ZfBZ  =  -46.9  +  20  log  Rkm .  (Back) 


D.  PRECIPITATION  PERFORMANCE 


Equation  (1)  can  be  used  to  compute  the  expected  value  of  the  received  signal  power  at 
the  antenna  output  port  for  reflectivities  of  different  levels  and  ranges.  Substituting  the 
radar  system  parameters  in  Eq.  (1),  one  obtains 

^  =  -80.1- 20  log  Rkm  +  ZdBZ  dBm.  (5) 

A  graph  of  received  weather  signal  power  vs.  range  for  reflectivities  between  5  and  75 
dBZ  is  shown  in  Figure  7.  The  curves  show  that  light  precipitation  will  produce  a  return 
at  or  above  the  system  noise  floor  for  ranges  to  45  km.  For  ranges  less  than  15-20  km., 
high  reflectivities  can  result  in  a  return  that  will  drive  the  receiver  into  saturation.  This 
situation  would  most  probably  occur  if  the  radar  were  used  to  observe  a  tornado  at  short 
range  as  funnels  are  known  to  exhibit  high  reflectivity.  It  should  also  be  remembered  that 
individual  pulses  can  exceed  the  average  power  level  by  as  much  as  10  dB  (See 
Appendix  1).  Thus,  saturation  is  an  issue  for  any  reflectivity  that  produces  average 
power  greater  than  a  threshold  set  10  dB  below  the  receiver  saturation  level  or  about  -40 
dBm.  (Back  to  Table  of  Contents) 
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Figure  7.  MWR-05XP  single  pulse  received  average  power  vs.  range  for  reflectivity  in 
the  range  5  to  75  dBZ.  Dashed  curves  show  the  system  noise  floor  and  receiver 
saturation  levels.  (Back) 


E.  BIRD  AND  INSECT  RETURN 


During  periods  of  clear  weather  or  in  areas  devoid  of  precipitation,  radar  echos  generally 
occur  due  to  clear  air  turbulence,  birds  or  insects.  It  is  not  always  easy  to  distinguish 
between  these  sources  of  return  [6].  However,  such  return  can  contaminate  data, 
especially  by  biasing  wind  velocity  measurements. 


The  basic  radar  range  equation  can  be  used  to  calculate  the  echo  power  from  insects.  For 
the  MWR-05XP,  the  signal  power  at  the  antenna  output  (receiver  input)  is  given  by 


Pr  =-74.7  +  10  log 


dBm. 


(6) 


The  scattering  characteristics  of  a  number  of  different  birds  and  insects  can  be  found  in 
[5].  Typical  X-band  radar  cross-sections  are  shown  in  Table  1. 
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Insect 

X-Band  RCS  in  cm2 

Aphid 

10'5 

Mosquito 

5x1  O'4 

Ladybug 

10'2 

Moth 

1.0 

Sparrow 

1.6 

Grackle 

16 

Table  1.  X-band  radar  cross-sections  for  some  birds  and  insects.  After  [5].  (Back) 


Using  the  RCS  values  in  Table  1,  the  echo  power  for  a  single  bird  or  insect  can  be 
calculated  for  various  ranges  using  Eq.  (6).  The  result  is  shown  in  Figure  8.  The  range 
of  RCS  in  Figure  8  is  10'  cm"  to  10  cm";  the  smallest  RCS  value  corresponding  to  an 
aphid  and  the  largest  to  a  big  bird.  It  is  clear  that  the  single  pulse  return  from  a  single 
bird  or  insect  will  generally  be  below  the  MWR-05XP  system  noise  level  and  is  not 


Figure  8.  MWR-05XP  single  pulse  echo  power  return  for  single  birds  and  insects.  (Back) 


a  concern.  However,  flocks  of  birds  and  insect  swarms  produce  a  much  higher  return  and 
such  returns  are  frequently  observed  during  the  summer,  both  day  and  night. 

The  return  for  flocks  of  birds  and  insect  swarms  is  best  characterized  in  terms  of  RCS 
density,  rj,  in  m  An  .  The  radar  range  equation  for  this  case  has  the  fonn 
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(V) 


PtfX^ejpJn^cv) 

r  1024(ln2)i?2  T1' 

Substituting  the  MWR-05XP  system  parameters  yields 

P  =  35  +  1 0  log  7(m2/cm3  -  20  log  Rkm  dBm.  (8) 

Typical  RCS  densities  for  birds  and  insects  have  been  summarized  by  Skolnik  [5].  Table 
2  shows  typical  values  of  RCS  density  at  X-band. 


Type  of  Flock  or  Swarm 

RCS  Density,  77Cm2/cm3 

Crows,  gulls,  geese,  ducks 

101U  -  10  s 

Blackbird  roost 

OS 

o 

l 

o 

Medium  sized  butterflies 

p 

c 

1 

O 

oc 

Aphid  migration  (major) 

io y-  i(T8 

All  insects  (1  hour) 

10  ^  —  10'6 

Table  2.  X-band  RCS  densities  for  some  bird  flocks  and  insect  swarms.  After  [5], 


(Back) 


MWR-05XP  Single  Pulse  Bird  Flock  or  Insect  Swarm  Return: 
RX  Power  in  dBm  vs.  Range  in  Kilometers 
RCS  Density  -100  dBcm2/cm3  to  -60  dBcm2/cm3 

Prof.  Jeffrey  B.  Knorr 

Naval  Postgraduate  School  Monterey,  CA  93943 


-  -  -  RX  Noise 
- •1.00E-10 

1.00E-08 

1.00E-06 
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Figure  9.  MWR-05XP  single  pulse  echo  signal  average  power  vs.  range  for  flocks  of 
birds  and  swarms  of  insects  with  various  RCS  densities.  Typical  RCS  densities  are  in  the 
range  10'10-  10"6  cm2/cm3.  (Back) 
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Figure  9  shows  the  result  obtained  when  the  RCS  densities  from  Table  2  are  substituted 
in  Eq.  (8).  The  curve  for  r)  =  10  6  cmVcnT  results  in  the  highest  return  that  might  be 
expected  and  corresponds  to  an  hourly  average  for  all  insects.  Birds,  although  they  have 
much  larger  individual  cross-sections,  have  much  lower  number  densities  and  therefore 
would  not  be  expected  to  produce  a  return  much  greater  than  that  shown  by  the  lowest 
curve  in  Figure  9,  77  =  10~10  cm2/cm3. 

The  conclusion  that  can  be  drawn  from  Figure  9  is  that  insect  swarms  can  produce  a 
significant  return  at  X-band.  Comparing  the  results  in  Figure  9  with  those  in  Figure  7 
leads  to  the  conclusion  that  insect  swarms  or  flocks  of  birds  can  produce  a  return  that  will 
be  displayed  as  a  reflectivity  in  the  range  15-55  dBZ.  WSR-88D  imagery,  available  on 
the  internet,  frequently  shows  clear  air  return  in  the  lower  part  of  the  range  predicted 
above.  It  should  also  be  noted  that  ornithologists  and  entomologists  use  radar  as  a  tool  in 
their  studies  of  birds  and  insects.  (Back  to  Table  of  Contents) 

F.  CLEAR  AIR  PERFORMANCE 

Scattering  from  clear  air  turbulence  can  also  be  observed  with  radar.  The  RCS  density 
for  clear  air  turbulence  is  [2] 


= o.39r»c;  (9) 

where  Cn  =  refractive  index  structure  parameter  in  nf  .  Typical  values  of  the  structure 
parameter  are C(2  «  6x10  1 7  (weak)  -  C2  «  3x  1 0  13  (strong)  [2].  Substituting  Eq.  (9)  in  the 
radar  equation  yields 


_  PtG;^ 
P  = 


>71 


;)(cr) 


1024(ln(2))R2 


(0.39)  C;. 


(10) 


Substitution  of  the  MWR-05XP  radar  system  parameters  then  yields  the  expression  for 
the  single  pulse  received  signal  power  due  to  scattering  from  clear  air  turbulence, 

^  =  15.9  +  101og(C,2)-201og(^„)  dBm.  (11) 


The  received  signal  power  for  strong,  medium  and  weak  turbulence  is  shown  in  Figure 
10.  The  figure  shows  that  the  return,  even  from  strong  from  clear  air  turbulence,  is  more 
than  20  dB  below  the  MWR-05XP  system  noise  floor.  This  return  is  much  weaker  than 
that  from  insects  or  precipitation.  Thus,  during  periods  when  no  precipitation  is  present, 
return  from  insects  and  birds  will  normally  be  much  greater  that  that  from  clear  air 
turbulence.  A  long  integration  time  would  be  required  to  enable  observation  of  clear  air 
turbulence  with  the  MWR-05XP.  Radars  designed  to  detect  clear  air  turbulence  normally 
operate  at  much  longer  wavelengths,  typically  10-20  times  greater. 
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MWR-05XP  Single  Pulse  Clear  Air  Performance: 
Received  Power  in  dBm  vs.  Range  in  Kilometers 
For  Turbulance  Parameter  CnA2  Strong,  Medium  &  Weak 
Prof.  Jeffrey  B.  Knorr 

Naval  Postgraduate  School  Monterey,  CA  93943 

0 

-20 
m  '40 

c  -60 
5  -80 

o  -100 

Q. 

T3  -120 
o 

|  -140 
-160 
-180 
-200 

Range  in  Kilometers 


Figure  10.  MWR-05XP  single  pulse  received  echo  signal  power  vs.  range  for  strong, 
medium  and  weak  clear  air  turbulence.  (Back) 

G.  FLUCTUATING  TARGET  INTEGRATION  GAIN 


-  -  -  RX  Noise 

Strong  IE-13 
Med.  IE-15 
Weak  IE-17 


Previous  figures  show  MWR-05XP  average  echo  power  for  a  single  pulse;  for 
precipitation,  birds,  insects  and  clear  air  turbulence.  Albersheim  has  published  closed 
form  expressions  relating  detection  probability,  Pd,  false  alann  probability,  Pfa,  pulses 
integrated,  n  and  required  SNR  for  non- fluctuating  targets  [7].  The  Albersheim’s 
equations  can  also  be  used  to  approximate  the  gain  for  incoherent  integration  of  n 
samples  [8].  Albersheim’s  expressions  are 


■4  =  ln[0.62/pJ 
5  =  In  [/>/(! -.(>)] 


(SNR\  =  xx  =-51og10« 

~  (  4.54 

6.2  + 


\ 


1CT7  <  Pfa  <  10  3 


v  y/n  +  0.44  j 

Q.\<Pd  <0.9, 


log10(X  + 0.12^5 +  1.7)  dB 
1  <  n  <  8096. 


(12) 


For  the  range  of  parameters  given  by  Albersheim,  the  error  in  the  estimate  of  (SNR)  i  is 
claimed  to  be  less  that  0.2  dB. 
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From  Eq.  (12),  it  follows  that  the  gain  for  incoherent  integration  of  n  pulses  is  given  by 

[8] 


k  =  A  +  0A2AB  +  1.1B 
f  0.454  A 


/(«)=  -0.38 

V  7  \y/n  +  0.44 ) 

G*c(n)  =  '%f(n)  0r 

Gnc{n)  =  5logl0n-  f  (n)l0logl0  k  dB. 


(13) 


The  MWR-05XP  noncoherent  integration  gain  for  a  detection  probability  of  0.95  and  a 
false  alarm  time  of  1  hour  is  shown  in  Figure  1 1 . 


Figure  1 1 .  MWR-05XP  postdetection  integration  gain  vs.  number  of  pulses  integrated 
for  Pd  =  0.95,  Tfa  =  1  hour  and  an  assumed  nonfluctuating  target.  (Back) 


The  primary  target  of  interest  is  precipitation,  which  produces  a  fluctuating  echo  signal 
return.  Thus,  further  analysis  is  required  to  determine  the  integration  improvement  for 
precipitation  echo  pulses. 
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Marshall  and  Hitschfeld  described  the  nature  of  the  fluctuating  return  from  precipitation 
[9].  The  number  and  diameter  of  drops  within  the  observation  volume  determine  the 
average  power  return.  The  return  power  for  a  particular  pulse  can  only  be  described 
statistically,  however.  Precipitation  is  a  Rayleigh  target  and  the  echo  signal  power  is 
exponentially  distributed  while  the  envelope  of  the  detected  signal  voltage  is  Rayleigh 
distributed.  Precipitation  is  thus  a  Swerling  target  of  Type  1  or  Type  2,  depending  on 
signal  correlation  [10].  Kanter  [11]  has  determined  the  detection  probability  for  a 
Rayleigh  target  when  pulses  are  arbitrarily  correlated  and  Shnidman  [12]  has  provided 
closed  form  equations  analogous  to  Albersheim’s  but  applicable  to  Swerling  targets  of  all 
classes.  Shidman’s  equations  are 


oo 

1 


K  =  \ 


n 


2 


[2  n 


Swerling  0 
Swerling  1 
Swerling  2 
Swerling  3 
Swerling  4 


JO  n  <  40 
[0.25  ft  >  40 


,  =  ^-0.81n(4^(l-^))+sgn(^-0.5)^-0.8ln(4i>(l-i>)) 


0=7 


Cj  =(l/if){[(17.7006P]  -18.4496)/5,  +14.5339]i>  -3.525} 

C2  =  (l/K)  {exp  (27.3  \Pd  -  25. 14)  +  (Pd  -  0.8)  [o.7  In (lO ~5/Pfa )  +  (2 n  -  20) /80]J 

_JC,  0.1  <Pd  <0.872 
dB~\c2+C2  0.872  <  Pd  <  0.99 

Q  —  ^q(Qb/10) 


SNR  =  Zl=(CXjn) 

SNRdB=  101og10(^) 

0.1  <7^  <0.99,  10'9  <  ^  <  10  3  1  <  ft  <  100. 


(14) 


These  equations  are  valid  for  a  larger  range  of  Pd  and  Pfa  but  a  smaller  range  for  n. 
Practically,  however,  once  n  exceeds  the  range  for  which  the  equations  are  valid, 
fluctuation  loss  is  negligible  and  integration  gain  can  be  determined  from  Eq.  (13). 

It  is  of  interest  to  investigate  the  single  pulse  fluctuation  loss  for  precipitation,  Swerling 
Case  1  (p=  1,  scan-to-scan  fluctuations)  or  Case  2  (p=  0,  pulse-to-pulse  fluctuations). 
This  can  be  determined  from  Shnidman’s  equations  and  is  illustrated  in  Figure  12. 
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Figure  12.  Fluctuation  loss  vs.  detection  probability  for  Swerling  Type  1  and  Type  2 
targets.  The  false  alarm  time  of  1  hour  corresponds  to  Pfa  =  6.3x10  10  for  the  MWR- 
05XP.  (Back) 


Fluctuation  loss  is  the  reduction  in  postdetection  integration  gain  relative  to  a 
nonfluctuating  target  if  only  a  single  pulse  is  available  to  make  the  detection  decision. 
Equivalently,  fluctuation  loss  is  the  increase  in  single  pulse  SNR  (relative  to  a 
nonfluctuating  target)  required  to  achieve  a  specified  detection  probability  when  the 
target  echo  signal  power  is  fluctuating  according  to  one  of  the  Swerling  models.  Figure 
12  shows  this  loss  is  12  dB  for  a  detection  probability  Pd  =  0.95  and  17.9  dB  for  a 
detection  probability  Pd  =  0.99. 

It  is  interesting  to  note  that  Figure  12  shows  the  single  pulse  fluctuation  loss  increases  as 
detection  probability  increases.  For  a  nonfluctuating  target  and  for  fixed  Pfa,  increasing 
single  pulse  detection  probability  also  requires  an  increase  in  SNR.  Since  the  single 
pulse  SNR  for  precipitation  (a  Rayleigh  target)  is  not  constant,  but  rather  is  exponentially 
distributed,  the  conclusion  to  be  drawn  is  a  larger  single  pulse  SNR  (relative  to  a 
nonfluctuating  target)  is  required  to  assure  a  single  precipitation  echo  is  detected  with 
high  probability.  The  opposite  argument  holds  for  small  detection  probability.  For  small 
detection  probability,  a  smaller  single  pulse  SNR  is  required  because  there  is  a  reasonable 
probability  that  a  large  echo  signal  return  from  precipitation  will  occur. 
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Fluctuation  loss  for  Swerling  Type  1  {p  =  1)  and  Type  2  (p  =  0)  targets  can  be  found 
from  Figure  12  or  from  Shnidman’s  equations.  According  to  Kanter  [1 1],  the 
improvement,  /  (/?) ,  for  incoherent  integration  of  n  pulses  with  correlation  coefficient  p 
is  given  by 


ne  =  min  [n,  1  +  (n  - 1)  In  (1/  /?)] 

/  ( n)  =  — — 

,U  (Lff'] 


(15) 


where  Lf  is  the  single  pulse  fluctuation  loss,  ne  is  the  equivalent  number  of  independent 
pulses  and  Gnc  is  the  gain  for  integration  of  n  uncorrelated  pulses,  Eq.  (13). 

It  is  of  interest  to  examine  the  effect  of  correlation  on  njn.  This  is  shown  in  Figure  13. 


ne/n  vs.  Correlation  Coefficient 
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Figure  13.  Effect  of  weather  signal  correlation  on  normalized  effective  number  of 
independent  pulses  for  one  to  twenty  pulses  received.  (Back) 


Figure  13  clearly  shows  that  when  the  echo  signal  sample  correlation  coefficient  exceeds 
p  ~  0.4 ,  the  effective  number  of  independent  pulses,  ne,  will  be  less  than  the  actual 
number  of  received  pulses,  n.  In  the  limit,  as  p  >  I  ,  the  effective  number  of 
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independent  pulses  converges  to  ne  =  1 ,  regardless  of  the  actual  number  received.  Thus, 

/  \(l/«e) 

for  a  large  number  of  received  pulses,  the  effective  fluctuation  loss,  Lfe  =  ( Lf  J  ,  will 
vary  from  its  single  pulse  value  when  p  =  1 ,  to  Lfe  » 1  (0  dB)  when  p  =  0  . 

The  overall  integration  improvement  for  a  fluctuating  target  can  be  determined  from  Eq. 
(15)  in  terms  of  the  number  of  pulses  and  their  correlation.  This  is  shown  in  Figures  14a 
-  14e. 
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Figure  14a.  Integration  gain  vs.  pulses  integrated  for  a  detection  probability,  Pd  =  0.95, 
false  alarm  probability  Pfa  =  e  .3xlO"10.  Pulse-to-pulse  correlation  is  p=  0.  (Back) 


A  conclusion  that  can  be  drawn  from  Figures  14a  -  14e  is  for  a  specified  correlation, 
p  <  1 ,  the  integration  improvement  approaches  that  for  a  nonfluctuating  target  (Swerling 

Case  0)  if  enough  pulses  are  integrated.  Equivalently,  ne  — >  n  and  the  fluctuation  loss, 

(\  \jne 

Lf  1  — >- 1  if  a  large  number  of  pulses  are  integrated.  The  convergence  occurs  quickly 

for  p  <  0.4  and  slows  as  p  — »  1 .  The  correlation  coefficient  depends  on  the  decorrelation 
time  for  the  weather  signal  and  the  radar  sample  interval.  Signal  correlation  is  discussed 
in  Ch.  III.  (Back  to  Table  of  Contents) 
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Integration  Gain  vs.  Pulses  Integrated 
With  and  Without  Fluctuation  Loss  (Rho  =  0.4) 
Pd  =  0.95,  Tfa  =  1  Hr. 
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Figure  14b.  Integration  gain  vs.  pulses  integrated  for  a  detection  probability,  Pd  =  0.95, 
false  alarm  probability  Pfa  =  6  .3xlO'10.  Pulse-to-pulse  correlation  is  p=  0.4.  (Back) 


Integration  Gain  vs.  Pulses  Integrated 
With  and  Without  Fluctuation  Loss  (Rho  =  0.6) 
Pd  =  0.95,  Tfa  =  1  Hr. 

Prof.  Jeffrey  B.  Knorr 

Naval  Postgraduate  School  Monterey,  CA  93943 


Figure  14c.  Integration  gain  vs.  pulses  integrated  for  a  detection  probability,  Pd  =  0.95, 
false  alarm  probability  Pfa=  6  .3xl0'10.  Pulse-to-pulse  correlation  is  p  =  0.6.  (Back) 
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Integration  Gain  vs.  Pulses  Integrated 
With  and  Without  Fluctuation  Loss  (Rho  =  0.8) 
Pd  =  0.95,  Tfa  =  1  Hr. 
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Figure  14d.  Integration  gain  vs.  pulses  integrated  for  a  detection  probability,  Pd  =  0.95, 
false  alarm  probability  Pfa=  6  .3xlO'10.  Pulse-to-pulse  correlation  is  p=  0.8.  (Back) 


Integration  Gain  vs.  Pulses  Integrated 
With  and  Without  Fluctuation  Loss  (Rho  =  0.99) 
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Figure  14e.  Integration  gain  vs.  pulses  integrated  for  a  detection  probability,  Pd  =  0.95, 
false  alarm  probability  Pfa=  6  .3xl0'10.  Pulse-to-pulse  correlation  is  p=  0.99.  (Back) 
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III.  TEMPORAL  AND  SPATIAL  SAMPLING  STRATEGIES 

A.  SIGNAL  CORRELATION 


A  weather  signal  is  typically  assumed  to  have  a  Gaussian  power  spectrum  with  spectral 
width  <jf  determined  by  the  rms  velocity  spread  crv  of  the  weather  system.  Spectral 
spread  and  velocity  spread  are  related  by  <j  f  =  2 cr,  /  A .  The  weather  signal 


autocorrelation  function  is  found  from  the  Fourier  transform  of  the  power  spectrum.  The 
correlation  coefficient  is  given  accordingly  as  (Refs.  [1]  -  [3]) 


P(T) 


(16) 


where 

t  =  time  in  seconds 

cr,  =rms  velocity  spread  in  m/s 


A  =. wavelength  in  meters. 


Figure  15  shows  the  expected  X-band  weather  signal  correlation  coefficient  for  parameter 
(Tv- 


Figure  15.  Correlation  coefficient  vs.  time  for  a  weather  signal  observed  at  f  =  9.7  GHz. 
RMS  velocity  spread  parameter  0. lm/s<  cr,  <10m/s.  (Back) 
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Of  interest  is  the  time  the  signal  remains  correlated  and  the  time  required  for  the  signal  to 
decorrelate.  Both  times  depend  on  RMS  velocity  spread,  cr,.  Velocity  spread  is  a 
function  of  turbulence  and  wind  shear. 

For  observations  in  X-band,  Figure  15  shows  the  signal  decorrelates  in  a  time  on  the 
order  of  5  -  10  milliseconds  when  RMS  velocity  spread  is  greater  than  1  m/s.  If 
p{j)  =  e  4  ~  -02  is  chosen  as  the  threshold  for  decorrelation,  then  decorrelation  time  can 
be  found  from  Eq.  (16)  as 

t  =  (/l/V2;rcrv).  (17) 

Figure  16  shows  the  decorrelation  time  obtained  by  applying  this  threshold  to  an  X-band 
radar  operating  at  9.7  GHz. 


Decorrelation  Time  vs.  RMS  Velocity  Spread 
for  a  9.7  GHz  Weather  Signal 
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Figure  16.  Decorrelation  time  vs.  RMS  velocity  spread  for  a  weather  signal  observed 
with  an  X-band  radar.  Threshold  for  decorrelation  is  assumed  /?( r)  =  exp (-4)  «  0.02  . 
(Back) 


Figure  16  shows  when  a  weather  target  has  RMS  velocity  spread  cr  ,  >  1  m/s ,  independent 
weather  signal  samples  can  be  obtained  with  an  X-band  radar  if  Ts  >  7  ms  where  Ts  is  the 
sample  time. 
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Also  of  interest  is  the  maximum  time  that  can  elapse  between  samples  if  it  is  desired  that 
the  samples  be  highly  correlated.  This  applies  to  velocity  estimation  where  the  pulse  pair 
algorithm  is  used.  If  y>(r)  =  e1'2  «  0.61  is  chosen  as  the  threshold,  then  signal  correlation 
time  can  be  found  from  Eq.  (16)  as 

r  =  /t/4;rcrv  seconds.  (18) 

Figure  17  shows  correlation  time  for  a  weather  signal  observed  with  an  X-band  radar. 


Correlation  Time  vs.  RMS  Velocity  Spread 
for  a  9.7  GHz  Weather  Signal 
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Figure  17.  Correlation  time  vs.  RMS  velocity  spread  for  a  weather  signal  observed  with 
an  X-band  radar.  Threshold  for  correlation  is  assumed  /?(r)  =  e~1/2  «  0.61.  (Back) 

Research  has  shown  that  90%  of  severe  storms  have  RMS  velocity  spread  cr,  <  8  m/s. 
Figure  17  shows  that  weather  signal  samples  separated  by  no  more  than  0.25  ms  will 
have  high  correlation  if  cr,  <  10  m/s .  Equivalently,  a  pulse  repetition  frequency, 

f  >  2.5  kHz  will  suffice  for  velocity  estimation  under  most  conditions.  At  the  same 
time,  Figure  16  shows  that  for  cr,  >  1  m/s  ,  7  ms  should  elapse  between  pulse  pairs  to 
assure  that  independent  sample  pairs  are  obtained.  An  even  longer  time  would  be 
desirable  to  accommodate  conditions  where  cr ,  <  1  m/s .  (Back  to  Table  of  Contents) 
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B.  SAMPLING  STRATEGY 


The  results  of  the  previous  section  are  the  basis  for  developing  an  optimum  sampling 
strategy.  Here,  “optimum  sampling  strategy”  means  the  strategy  that  will  both  maximize 
volumetric  update  rate  and  minimize  the  variance  of  weather  signal  parameter  estimates. 
Equivalently,  it  is  the  strategy  that  will  result  in  the  most  rapid  collection  of  independent 
sample  pairs  for  each  of  the  many  radar  beam  positions  (and  resolution  cells)  of  interest. 

Sampling  and  volumetric  update  requirements  appear  conflicting.  A  pair  of  closely 
spaced  pulses  is  required  to  estimate  velocity  (and  avoid  velocity  ambiguity).  As  seen 
above,  spacing  should  not  exceed  0.25  ms.  in  X-band.  Close  spacing  of  pulses  is 
consistent  with  a  rapid  volumetric  update  rate.  However,  a  large  number  of  independent 
sample  pairs  is  required  to  estimate  average  velocity  with  small  variance  and  this  means 
successive  sample  pairs  should  be  widely  spaced  in  time.  Successive  sample  pairs  should 
be  separated  by  at  least  7  ms.  in  X-band;  longer  if  velocity  spread  is  less  than  1  m/s.  For 
a  mechanically  scanned  radar,  these  requirements  indeed  result  in  a  conflict.  For  an 
electronically  scanned  radar,  they  do  not.  Rapid  volumetric  update  and  parameter 
estimates  with  small  variance  can  be  simultaneously  achieved. 

The  MWR-05XP  uses  a  combination  of  frequency  and  phase  control  to  steer  the  radar 
antenna  beam  electronically.  Switching  time  for  the  ferrite  phase  shifters  is  in  the  range 
90  -  100  ps.  Oscillator  settling  time  is  shorter.  Thus,  beam  steering  is  constrained  by 
phase  shifter  switching  time.  It  will  be  assumed  here  that  100  ps  is  required  to  change 
beam  positions. 

The  MWR-05XP  pulse  repetition  frequency  is  limited  by  the  average  output  power 
capability  of  the  TWT  used  in  the  transmitter.  On  a  continuous  basis,  the  maximum  duty 
cycle  of  approximately  0.01  limits  PRF  to  10  kHz.  This  PRF  is  high  enough  to  achieve  a 
pulse-to-pulse  correlation,  p  ( 7) )  >  0.61 ,  while  accommodating  velocity  spreads  as  small 

as  crv  =  0.25  m/s  in  X-band.  The  pulse  repetition  interval  for  this  PRF  is  100  ps,  the  same 
as  the  switching  time  for  the  phase  shifters.  Thus,  the  time  required  to  transmit/receive 
two  pulses  and  switch  beam  positions  is  300  ps. 

Consider  a  raster  scan  consisting  of  200  beam  positions.  If  a  single  pair  of  pulses 
separated  by  0.25  ms.  was  transmitted  in  each  of  the  200  beam  positions,  an  elapsed  time 
of  60  ms  would  be  required  to  scan  the  frame.  Figure  16  shows  that  repeating  a  scan  of 
the  frame  after  an  elapsed  time  of  60  ms.  would  result  in  an  independent  sample  pair  even 
for  a  velocity  spread  as  small  as  cr,  =0.1  m/s  .  This  rapid  frame  scan  time  enables 

collection  of  a  sufficient  number  of  independent  samples  for  computation  of  parameter 
estimates  with  small  variance  on  a  time  scale  consistent  with  requirements  for 
investigation  of  rapidly  developing  storm  features.  (Back  to  Table  of  Contents) 

C.  DATA  QUALITY 


The  sampling  strategy  described  above  allows  for  simple  control  of  data  quality.  The 
desired  variance  of  a  parameter  estimate  can  be  specified  and  since  samples  are 
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independent,  the  number  of  samples  required  can  be  easily  determined.  Collecting  the 
computed  number  of  samples  with  sufficiently  high  SNR  will  then  result  in  weather 
signal  parameter  estimates  with  known  variance. 

For  example,  suppose  it  is  desired  to  estimate  reflectivity,  Z  ,  with  standard  deviation 
(, SD ),  oy  <  +0. 1Z  .  Z  is  estimated  from  echo  signal  average  power,  P  .  The  standard 

deviation  for  a  single  sample  of  the  signal  power  is  P  (See  Appendix  1  and  Ref.  [9]).  For 
n  independent  samples  the  standard  deviation  is  p/ 4n  .  n  =100  samples  will  reduce  the 

standard  deviation  of  the  average  power  estimate  to  SD  |pj  =  0.  IP ,  thereby  yielding  the 

desired  SD  of  the  reflectivity  estimate,  Z  .  In  dBZ,  this  100  sample  reflectivity  estimate 
has  a  standard  deviation  of  approximately  Sd[z j  «  ±0.46  dBZ  .  The2cr  accuracy  of  the 

estimate  is  better  than  ±1  dBZ  .  By  the  Central  Limit  Theorem,  the  probability  the  n  = 
100  sample  reflectivity  estimate  is  within  2crof  the  true  reflectivity  is  0.95. 

The  previous  example,  estimating  reflectivity  with  specified  variance,  requires  that  echo 
signal  power  be  accurately  measured.  Thus,  in  the  case  of  reflectivity,  the  radar  must  be 
accurately  calibrated  to  achieve  high  data  quality.  (Back  to  Table  of  Contents) 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

Electronic  scan  gives  the  MWR-05XP  radar  unique  capabilities  as  a  scientific  instrument 
for  investigation  of  weather  phenomena.  This  report  provides  performance  curves  that 
quantify  the  signal  return  expected  from  precipitation  with  varying  intensity  and  in  clear 
air,  the  return  expected  from  birds,  insects,  atmospheric  turbulence  and  wind  shear. 
Curves  are  also  provided  to  determine  postdetection  integration  gain  when  fluctuation 
loss  occurs  and  pulses  are  correlated.  Lastly,  sampling  requirements  are  discussed  in  the 
context  of  weather  signal  correlation.  It  is  shown  that  electronic  scan  enables  collection 
of  independent  samples  while  simultaneously  maintaining  a  rapid  volumetric  update  rate. 
This  is  accomplished  by  collecting  a  single  sample  pair  for  each  beam  position  in  a 
frame,  then  repeating  the  scan  of  the  frame  to  obtain  another  set  of  sample  pairs  that  will 
be  independent  of  the  first  set  due  to  the  signal  decorrelation  that  occurs  during  the  frame 
scan  time.  Both  the  sampling  scheme  and  the  use  of  frequency  for  beam  scanning  present 
new  signal  processing  challenges.  (Back  to  Table  of  Contents) 

B.  RECOMMENDATIONS 

There  are  several  areas  in  need  of  future  investigation. 

(1)  The  MWR-05XP  must  be  calibrated  to  obtain  accurate  estimates  of  reflectivity. 
The  various  approaches  to  calibration  need  to  be  investigated  and  at  least  one 
approach  needs  to  be  developed  and  implemented  for  use  with  the  MWR-05XP. 

(2)  The  optimum  sampling  strategy  described  in  this  report  results  in  data  files  with  a 
different  structure  than  those  produced  by  a  mechanically  scanned  radar  antenna. 
Software  must  be  developed  to  process  the  data  in  these  files  and  compute 
estimates  of  reflectivity,  average  velocity  and  velocity  spread. 

(3)  A  third  area  in  which  interesting  work  could  be  done  is  data  quality.  It  should  be 
possible  to  specify  a  desired  variance  and  confidence  level  for  signal  parameter 
estimates  and  use  these  to  drive  the  sample  collection  process.  This  would  lead  to 
estimates  of  reflectivity,  average  velocity  and  velocity  spread  with  quantified 
variance. 

(Back  to  Table  of  Contents) 
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APPENDIX  1.  WEATHER  SIGNAL  PDF 


The  fluctuating  echo  signal  from  precipitation  has  been  studied  by  Marshall  and 
Hitschfeld  [9],  The  pdf  of  the  weather  signal  radar  cross-section  and  thus  the  echo 
signal  power  is  exponential, 


p(P)  =  =  e 


(A-l) 


The  pdf  of  log(P)  is  given  by 


where 


p{\ogP) 


exp 


f  log  P  exp  (log/3)} 
M  P  J 


MP 


M  =  log10  e  =  0.434  .  (A -2) 


Eq.  (A-2)  is  the  basis  for  a  graph  of  the  pdf  of  log10  ( P / Pj ,  shown  in  Figure  18.  Of 

interest,  is  the  range  over  which  the  precipitation  echo  signal  power  can  be  expected 
to  vary.  From  Figure  18,  it  is  clear  that  the  instantaneous  power  will  typically  lie  in 
the  range  from  30  dB  below  the  average  to  10  dB  above  the  average.  From  Eq.  (A-l), 
the  probability  the  power  exceeds  the  average  by  6  dB  is  only  0.018  or  slightly  less 
than  1.8%.  The  probability  the  power  exceeds  the  average  by  10  dB  is  almost 
negligible.  Thus,  if  a  threshold  is  set  10  dB  below  the  receiver  saturation  level, 
precipitation  return  with  average  echo  signal  power  at  or  below  this  threshold  will  not 
produce  any  echos  that  drive  the  receiver  into  saturation. 


For  the  MWR-05XP,  the  threshold  described  above  would  be  about  -40  dBm.  From 
Figure  7,  it  is  evident  that  if  reflectivity  is  below  55  dBZ,  the  receiver  will  operate 
linearly  for  all  ranges  greater  than  5  km.  (Back  to  Table  of  Contents) 
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PDF  of  Precipitation  Echo  Signal 

Prof.  Jeffrey  B.  Knorr 

Naval  Postgraduate  School  Monterey,  CA  93943 


(PIP avg)  in  dB 


Figure  18.  PDF  of  precipitation  echo  signal  power.  P  is  the  instantaneous  power  and 
Pavg  is  the  expected  value  of  the  echo  signal  power.  (Back) 
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APPENDIX  2.  MWR-05XP  SYSTEM  PARAMETERS 


Effective  Radiated  Power  (measured) 

Pulse  Width  (measured) 

PRF 

Frequency 
Azimuth  Beamwidth 
Elevation  Beamwidth 
Phase  Shifter  Switching  Time 
System  Noise  Floor  (measured) 

Receiver  Saturation  Level  (measured  at  RX  input) 

(Back  to  Table  of  Contents) 


109.7  dBm 
1.0  ps 

10  kHz  (max.) 
X-band 
1.8  degrees 
2.0  degrees 
100  ps 
-1 12.8  dBm 
-30  dBm 
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